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It has recently been suggested '^'^ that galaxies in the early Universe can grow through the 
accretion of cold gas, and that this may have been the main driver of star formation and 
stellar mass growth^ SIS Because the cold gas is essentially primordial, it has a very low 
abundance of elements heavier than helium (metallicity). As it is funneled to the centre of a 
galaxy, it will lead the central gas having an overall lower metallicity than gas further from 
the centre, because the gas further out has been enriched by supernovae and stellar winds, 
and not diluted by the primordial gas. Here we report chemical abundances across three 
rotationally-supported star-forming galaxies at z~3, only 2 Gyr after the Big Bang. We find 
an 'inverse' gradient, with the central, star forming regions having a lower metallicity than 
less active ones, opposite to what is seen in local galaxies l^l^. We conclude that the central gas 
has been diluted by the accretion of primordial gas, as predicted by 'cold flow' models. 
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We selected three Lyman-break galaxies^ among the AMAZE and LSD^'^sample, SSA22a- 
C16 (z = 3.065), CDFa-C9 (z = 3.219) and SSA22a-M38 (z = 3.288), to test the presence 
of metallicity gradients in galaxies beyond z = 3, when the Universe was only 14-16% of its 
actual age, and only a small fraction (~ 5 — 10%'^ of the total stellar mass in the Universe had 
been assembled. This epoch is particularly interesting, as it precedes the peak of the cosmic star 
formation density We selected these three galaxies because they show a remarkably symmetric 
velocity field in the [OIII] emission line, which traces the ionized gas kinematics (see Fig. [T]). Such 
kinematics indicates that these are rotationally supported disks (Gnerucci et al., in preparation), 
with no evidence for more complex merger-induced dynamics. Assuming an exponential disk mass 
model to reproduce the dynamical properties of the three sources, we derive dynamical masses in 
the range ~ 10^° — 2 10^^ Mq. By using the H/3 line flux 02, we infer star formation rates as high 
as ~ 120 Mq/ut after correcting for dust extinction!^. We carefully checked for the presence of 
Active Galactic Nuclei (AGN) in our galaxies, as they would invalidate the calibrations derived 
for star-forming galaxies because of the different ionizing continua involved (see Supplementary 
information). 

Near-infrared spectroscopic observations of the galaxies were obtained with the integral field 
spectrometer SINFONI ^ on the Very Large Telescope of the European Southern Observatory. 
The spectra cover the 1.45 — 2 Alum wavelength range; strong rest frame optical emission lines of 
the observed galaxies, such as [OIII]A5007, [OIII]A4959, H/3, [NeIII]A3870 and [OII]A3727, are 
redshifted to this range. The flux ratio between the (rest-frame) optical nebular lines observed in 
SINFONI spectra depend on gas metallicity, and have been cahbrated for local low and high metal- 
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licity galaxies^. In particular, we use a combination of three independent metallicity diagnostics 
^ [OIII]A5007/H/3, [OIII]A5007/[OII]A3727 and, when available, [NeIII]A3870/[OII]A3727, to 
derive the gas phase metallicity in our spectra. Even if a single line ratio can constraint the metal- 
licity in a broad interval of ~ 0.4 dex or more, the use of three independent diagnostics helps us 
in removing the degeneracies and reducing the uncertainties; we can determine for each spatial 
pixel the combination of metallicity and extinction that best fit the observed line ratios with our 
calibrations (see Supplementary information for the details). 

The derived metallicity maps are shown in Fig. [T| An unresolved region with lower metallic- 
ity is evident in each map, surrounded by a more uniform disk of higher metal content. In CDFa-C9 
the lower metallicity region is coincident with the galaxy center, as traced by the continuum peak, 
while it is offset by ~ 0.60" (4.6 kpc) in SS22a-C16 and ~ 0.45" (3.4 kpc) in SS22a-M38. On 
the other hand, in all the galaxies the area of lower metallicity is coincident or closer than 0.25" 
(1.9 kpc, which is half the point spread function full- width at half-maximum) to the regions of 
enhanced line emission, which trace the more active star forming regions. The average difference 
between high and low metallicity regions in the three galaxies is 0.55 in units of 12-i-log(0/H), 
larger than the ~ 0.2 — 0.4 dex gradients measured in the Milky Way and other local spirals'^ on 
the same spatial scales. The measured gas phase abundance variations have a significance between 
98% and 99.8% (see Tab.[T]). 

The observed line ratios depend on metallicity both directly and indirectly through the ion- 
ization parameter U (the ratio between the flux of ionizing photons and the gas density), with lower 
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metallicity star forming regions generally have higher U'^. In principle, the observed variation of 
line ratios across galaxies could still be due to variations of U unrelated to the metallicity distri- 
bution. As local galaxies show a broad range of ionization parameter, this dependence on U is 
responsible for the spread included in our calibrations. In order to exclude the possibility that the 
line ratios in our galaxies are biased due to significantly different values of ionization parameter at 
high redshift, we have examined their dependence on U based on the latest photoionization models 
available These models are not able to simultaneously reproduce the different observed 
hne ratios with a variation of the ionization parameter alone: at constant metallicity the observed 
line ratios differ from models predictions by ~ 0.2 dex. This demonstrate that different metallic- 
ities are required to reproduce the observations, and even if other effects are present, they are not 
dominant. 

Observed metalhcity gradients in local disk galaxies show metallicities decreasing outward 
from the galactic centers. In contrast, our high redshift sources show "inverted" gradients, with 
lower metallicity in a well defined active region closer to the galaxy center. Current models of 
chemical enrichment in galaxies SHIIIII cannot reproduce our observations at the moment, as they 
assume radially isotropic gas accretion onto the disk and the instantaneous recycling approximation 
Nevertheless, the detected gradients can be explained in the framework of the cold gas accretion 
scenario!^ recently proposed to explain the properties of gas rich, rotationally supported galaxies 
observed at high redshift '2^'^. In this scenario, the observed low metallicity regions are created 
by the local accretion of metal-poor gas in clumpy streams ^; this gas penetrate deep into the 
galaxy following the potential well, and sustaining the observed high star formation rate in the pre- 
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enriched disk. Stream-driven turbulence is then responsible for the fragmentation of the disks into 
giant clumps, as observed at z > 2^^, that are the sites of efficient star formation and possibly 
the progenitors of the central spheroid. This scenario is also in agreement with the dynamical 
properties of our galaxies, which appear to be dominated by gas rotation in a disk with no evidence 
of the kinematical asymmetries typically induced by mergers . 

If infall and outflows are not present, metallicity is directly related to the fraction of gas 
present in the galaxy. For this reason additional hints on the presence of infall and outflows can 
be obtained considering these two quantities. We can derive an estimate of gas fractions in each 
single region of our galaxies assuming that the relation between SFR and gas density observed 
in the local universe (the Schmidt-Kennicutt relation also applies to high redshift galaxies 
The gas fractions obtained for the two regions in each galaxy are shown in Fig. |2] plotted againts 
metallicities Z. It is evident that our results do not support the closed box model (black line) 
where galaxies have no contribution from infalls and outflows. Our data can be reproduced only 
by including either outflows of very enriched material, such as the ejecta from SNe, which extract 
metals and decrease Z, or infalls of pristine gas which increase the gas fraction and decrease Z. In 
particular, if strong infall of pristine gas is present, we would expect higher gas fractions in the low 
metallicity regions, as observed. Using simple modelP^, we estimate the amount of infall and 
outflow required to reproduce the observations. We take into account two opposite scenarios. If 
inflows are dominant, we require ~ 3.5 times more pristine gas infall in the low metallicity regions 
than in the high metallicity part. On the other hand, if outflows are dominant, high metallicity 
regions can be explained with outflow ~ 2 times the SFR, while unrealistic outflows ~ 12 times the 
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SFR would be needed to explain the properties of low metallicity regions. We therefore conclude 
that the differences in gas fractions as derived by our estimates are not only reasonable given the 
measured metallicity differences, but are also consistent with expectations from massive infalls of 
low metallicity gas in the central regions of the galaxies. They are unlikely to be explained using 
outflows only. 

Our observations of low metallicity regions in three galaxies at z ~ 3 therefore provide the 
evidence for the actual presence of accretion of metal-poor gas in massive high-z galaxies, capable 
to sustain high star formation rates without frequent mergers of already evolved and enriched sub- 
units. This picture was already indirectly suggested by recent observational studies of gas rich 
disks at 2; ~ 1 — 2^ ^ by the study of the variation of metallicity as a function of gas fraction and 
star formation rate The data presented here are now directly revealing the metal poor accreted 
gas responsible for fueling the star formation rate in massive galaxies in the early Universe. 

Although the three galaxies discussed here are the only 2; ~ 3 sources for which a metallicity 
gradient has been measured at the moment, the detection of lower metallicity clouds in all three 
galaxies suggests that metal poor gas accretion is probably a common feature in actively rotation- 
aUy supported star forming galaxies at this early cosmic epoch. However, it is hard to speculate 
how much this accretion mechanism is dominant in the global star forming population, as the rel- 
ative number of rotationally supported and more dynamically complex galaxies is still unknown, 
and similar observations are still missing in different classes of galaxies. 
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Table 1 : Properties of the sample galaxies, z is the source redshift, the stellar mass derived from 
multi band photometry and spectral energy distribution fitting assuming a Chabrier'^ Initial Mass 
Function ^, while the last three columns show the metallicity integrated for the whole galaxy, for the 
low metallicity regions and for the high metallicity regions respectively, in units of 12+log(0/H). 
The amplitude and significance of the detected gradients were evaluated extracting ID spectra 
integrated over apertures comparable to the PSF (FWHM=0.5") in both the lowest and highest 
metallicity region, and using them to measure the line ratios (see Supplementary information). 
The uncertainties are the la confidence intervals. 
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Figure 1: Surface brightness and velocity of the [OIII]A5007 Hne, and metallicity maps. The data for SSA22a-C16 
(z = 3.065, upper panels), CDFa-C9 (z — 3.219, central panels) and SSA22a-M38 {z — 3.288, lower panels) were 
obtained with the SINFONI spectrograph using the 0.125 x 0.250" pixel scale in seeing-limited conditions, resulting 
in a spatial resolution of ^ 0.5" (Full Width Half Maximum of the Point Spread Function, PSF). The maps were 
extracted from the SINFONI datacube after a Gaussian smoothing with FWHM=3 pixels (0.375"). The left panels 
in each row show normalized surface brightness in the [OIII]A5007 emission line. The same line has been used to 
derive the velocity maps shown in the central panels; the observed gas kinematics is compatible with a rotating disk, 
with no evidence of merger-induced complex dynamics. The right panels show maps of gas phase metallicity, as 
relative abundances of oxygen and hydrogen parameterized in units of 12 + log{0 / H). Lower metallicity region, 
corresponding to a higher [OIII]/[OII] ratio, are surrounded by a more enriched disk. The crosses in each panel mark 
the position of the continuum peak. 
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Figure 2: Gas fractions for the different regions in the observed galaxies. To estimate gas masses in the low metal- 
licity regions and in the rest of the galaxies, we used the H/3 flux corrected for extinction. We convert measured star 
formation rate densities from H/3 to a gas mass by inverting the Kennicutt-Schmidt lawl^. To estimate stellar masses, 
we have to relay on the continuum flux from our SINFONI data, as we are still missing deep and high resolution 
continuum images of the sources. In the datacubes of our galaxies the continuum is detected, although at a lower 
S/N than the emission lines, making this estimate more uncertain. We have to further assume that the stellar mass 
scales as the K-band continuum light (V band rest frame). We thus measure the total K band continuum flux on the 
whole galaxy and the flux in the lower metallicity regions, using these relative flux measurements to rescale the total 
stellar mass from Spectral Energy Distribution fitting^ and estimate the stellar mass in the metal poor and metal rich 
regions. We plot the derived gas fraction against the metallicity in solar units Z/Zq. Filled blue and red circles show 
respectively the low and high metallicity regions of each galaxy; black dashed line, expectations from a closed-box 
model. The difference between the regions can be explained by a simple model including only infaU of gas fi, which 
is, respectively, 3.5 (blue solid line) and 1 times (red solid line) the SFR for the regions of lower and higher metal- 
licity. Dashed cyan lines, a model of pure outflow /o, 12 times the SFR. The observed difference in gas fraction it is 
consistent with the expectation from a massive infall of low metallicity gas in the central regions. The error bars are 
la formal uncertainties on the measuraments, given the assumptions discussed. 
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Supplementary Information 

Measuring metallicities: At the redshift of the galaxies in our sample (z = 3.0 — 3.3), five promi- 
nent rest- frame optical emission lines are observable in the H-i-K SINFONI grism: [OII]A3727 
and [NeIII]A3870 in the H band, H/3, [OIII]A4959 and [OIII]A5007 in the K band. The resulting 
line ratios allow us to measure the gas phase metallicity using five well calibrated diagnos- 
tics, i?23 = ([OII]A3727 -1- [OIII]A4959 -i- [OIII]A5007)/H/?E1, [OIII]A5007/H/3, [OII]A3727/H^, 
[Oin]A5007/[OII]A3727 and [NeIII]A3870/[OII]A3727. Although these calibrations were ob- 
tained in local galaxies, where the physical conditions, such as the ionization parameter, might 
be different respect to higher redshift'^, these strong line metallicity diagnostics do not deviate at 
high-z more than about 0.1 dex'^S. 

These calibrations only apply if the ionizing source is stellar radiation, while the presence 
of an AGN ionizing continuum would invalidate the calibrations derived for star- forming galaxies. 
This could be a problem for our sources, as the larger metallicity differences are detected close to 
the galaxy nucleus, where an eventual AGN contribution would be more important. We therefore 
checked that, for our galaxies, rest frame UV spectra, soft and hard X-ray data, as well as 24yum 
Spitzer data show no indication of AGN ^. Moreover, we verified with spectral classification 
diagrams'^ that the line ratios detected in the central nuclear regions of our sources are compatible 
with a pure starburst in a low metallicity galaxy, ruling out the presence of any AGN contribution 
to the detected gradients, or significantly different physical conditions in the star forming regions 
with respect to local galaxies. 
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Obviously, not all of the five diagnostics listed above are independent, and each has its 
own advantages and disadvantages. [OIII]A5007/H/3 is essentially unaffected by dust extinc- 
tion, but it allows two different metallicities for a given ratio (see Fig. [3]). On the other hand, 
[Oin]A5007/[On]A3727 has a monotonic dependence on metallicity, but it is affected by red- 
dening. [NeIII]A3870/[OII]A3727 is instead both monotonic and less affected by dust, but the 
[NeIII]A3870 line is much fainter than the others, and thus not always observable. However, the 
combination of these three indicators allows to measure unambiguously the metallicity, taking 
into account the effect of dust reddening (see Fig. |4]). We exclude the remaining two line ratios 
from the following analysis, as i?23 is very similar to [Oin]A5007/H/3 at low metallicities and 
[OII]A3727/H/3 is not independent from the others. 

For each spatial pixel of the datacube with S/N > 3 in [OIII]A5007/H/3, we determine the 
best combination of metallicity and extinction that minimizes the in the three selected diagnos- 
tics, after a Gaussian smoothing of the data with a FWHM=3 pixels (0.375", lower than the spatial 
resolution of ~ 0.5"). To compute the x^, we take into account both the uncertainties on the line 
ratio measurements and the dispersion of each calibration Such a dispersion is mostly due to the 
broad range of ionization parameters in the galaxies used to derive the calibrations^^. There is no 
indication that the metallicity calibrations and the spread based on local galaxies may fail at high 
redshift, because the measured line ratios for our galaxies are consistent with the ones measured 
in local galaxies, and the comparison with the photoionization models available do not support the 
presence of unusually high ionization parameters. The spread in the calibrations results in a large 
uncertainty when the metallicity is measured using a single line ratio (~ 0.2 dex or more, depend- 

15 



ing on the line ratio), but the combination of the three independent tracers helps us to decrease the 
formal uncertainty (see Fig. |5]) . The additional systematic uncertainty on the absolute metallicity 
associated with the photoionization models used to derived the calibrations is not considered, as 
we are only interested in relative metallicity differences. The fitting results are shown in Fig.|4]for 
the three galaxies. 

To quantify the significance of the detected gradients we extracted ID spectra integrated over 
circular apertures within a diameter of 0.5", comparable to the PSF. We selected two apertures in 
each galaxy, one centered on the lowest metallicity region and one with highest metallicity, and 
we used the resulting line ratios to derive the metallic ities in both regions. The fits are shown in 
Fig. [3] for the three galaxies, and the results are listed in Tab. 1. The allowed values for the dust 
extinction are in the range < Ay < 1.2, in agreement with the estimate derived from the SED 
fitting Ay ~ 0.7 — 1 for the three galaxies^. Although the extinction has a large uncertainty, the 
metallicity is relatively well constrained. 
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Figure 3: Diagnostic tools used to measure metallicities. The low and high metallicity regions for 
the three galaxies are plotted, SSA22a-C16 on the left, CDFa-C9 in the center, SSA22a-M38 on the 
right. The first panels show the [OIII]A5007/[OII]A3727 line ratio as a function of metallicity, the 
second panels show the [OIII]A5007/H/3, that is very similar to R23 (shown in the third panels but 
not used for the minimization) at low metallicities, where the [OII]A3727 flux is small compared 
to [OIII]. The fijurth panels show the [NeIII]A3870/[OII]A3727 ratio. The blue points are relative 
to the lower metallicity regions, and the red points to the higher metallicity ones. The errorbars 
denotes the 1 a uncertainties in the flux measurements or in metallicity fitting. An upper hmit is 
shown where the [Nelll] emission line was not detected. 
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Figure 4: Line ratios and metallicity maps for the three galaxies of the sample. These maps are 
obtained after a Gaussian smoothing with FWHM=3 pixels (0.375"). The left panels show the 
[Oin]A5007/[On]A3727 line ratio maps, while the [OIII]A5007/H/3 maps are presented in the 
central panels. The line ratios gradients have a significance between 3 and 5(t. The right panels 
show the resulting maps of gas phase metallicities, as relative abundances of oxygen and hydrogen 
parameterized in units of 12 + log{0/H). The crosses in each panel mark the position of the 
continuum peak. 
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Figure 5: Example of dependence on metallicity using three independent diagnostics. The 
variation of as a function of metallicity is shown as an example for the integrated spectra of 
SSA22a-M38 using the [OIII]/[OII] line ratio only (blue curve), [0II]/H/3 only (red curve), and 
[NeIII/[OII] only (magenta curve). The total from the combined three independent measure- 
ments is shown as a black curve: the global 1 a confidence interval, represented as a gray shaded 
region is significantly narrower in this case. The best fitting value is marked with a green dotted 
line. 
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